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Transient and steady state currents through dc-biased quantum impurity models beyond the linear 
response regime are of considerable interest, both from an experimental and a theoretical point of 
view. Here we present a new analytical approach for the calculation of such currents based on the 
flow equation method (method of infinitesimal unitary transformations). Specifically, we analyze 
the Anderson impurity model in its mixed valence regime where the coupling to the leads is switched 
on suddenly at time t = 0. We observe the real time buildup of the current until it reaches its steady 
state limit. 
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I. INTRODUCTION 

Transport properties of quantum devices beyond the 
linear response regime have generated a lot of interest 
in the past decade. Experimentally, this is due to the 
recent advances in nanotechnology that permit to ap- 
ply large electrical fields in low dimensional electronic 
structures. Theoretically, transport beyond the linear re- 
sponse regime is interesting since it explores genuine non- 
equilibrium quantum many-body phenomena. A particu- 
larly well-studied case, both experimentally and theoret- 
ically, are quantum dots in the Coulomb blockade regime 
that display Kondo physics [H-0| : here the shot noise 
generated by the steady state current serves as a source 
of decoherence that suppresses the Kondo quasiparticle 
resonance for sufficiently large voltage bias [4(, thereby 
reducing the differential conductance [5|. 

However, the interplay of correlation physics and non- 
equilibrium is difficult to address theoretically, in spite 
of considerable effort in recent years. New numerical 
methods have been developed like the scattering state nu- 
merical renormalization group [6] , Monte Carlo methods 
0,9; the time-dependent density matrix renormalization 
group and other real time methods [l3|, [H| • Ana- 

lytical approaches are perturbative Keldysh calculations 
jl5| , extensions of the renormalization group [3, fl6H28| , 
generalizations of NCA (non-crossing approximation) to 
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non-equilibrium |29H31f . l/7V-expansions [32], Gutzwillcr 
methods 1331 and various approaches builing on integra- 
bility [34H37]. Since all of these methods have their re- 
spective limited domain of applicability, there is still a 
need for new theoretical methods. 

In the past few years the flow equation method 
(method of infinitesimal unitary transformations) [38l.l39j 
was used for a number of non-equilibrium quantum 
many-body problems like interaction quenches [40l - |42l | 
and dc-transport beyond the linear response regime 



(la . [i7l . |43j . In particular, for the Kondo model numer- 
ous quantities like spin susceptiblity, magnetization and 
T-matrix have been calculated for large voltage bias in 
the steady state [13, S3] ■ In addition, the flow equation 
method is particularly suitable for calculating the real 
time evolution of non-equilibrium problems [44j . There- 
fore it offers the possibility to study the transient time- 
dependent buildup of a quantity until it reaches its steady 
state value, see for example the calculation of the mag- 
netization dynamics in the ferromagnetic Kondo model 
40]. This defines the question investigated in this pa- 
per: Using the flow equation method, we calculate the 
time-dependent buildup of the electrical current through 
an Anderson impurity model when the coupling between 
the leads and the quantum dot is suddenly switched on at 
time t = 0. Thereby we develop a new analytical method 
for calculating transport properties of interacting quan- 
tum systems beyond the linear response regime, both for 
transient and steady state behavior. 

The model of a single level quantum dot coupled to two 
leads is described by the Anderson impurity Hamiltonian: 
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k denotes the wave vector, a =f , I the electron spin and 
a — L,R the left and right lead. For time t < both 
leads are in equilibrium at different chemical potentials 
/jLl and hr. The hybridzation V between leads and the 
dot is then switched on at time t = and we are in- 
terested in the current I(t) as a function of time. For 
simplicity we restrict ourselves to symmetric coupling to 
the leads, although the calculation can be generalized in 
a straightforward way. 

An explicit expression for I(t) is achieved via the 
forward-backward technique of the flow equation method 
fjij : The current operator is expressed in the diagonal 
basis of the Hamiltonian (|T|), where its time evolution can 
be worked out easily. Then the time-evolved operator is 
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transformed back into the original basis, where the initial 
condition of non-interacting Fermi gases with different 
chemical potentials is given. This yields the final answer 
with an explicit expression for the current as a function 
of time. Approximations enter during the diagonalization 
step of the Hamiltonian, which limits our calculation to 
weak and intermediate interaction U. However, the volt- 
age bias can be large (beyond the linear response regime) 
and the real time evolution followed into the asymptotic 
steady state limit without any difficulties. 




FIG. 1: Schematic representation of the parameters in the 
Anderson impurity model. 



II. TRANSFORMATION OF THE 
HAMILTONIAN 

We employ a symmetric/antisymmetric basis Ck± a 
■^(ckLa ± CkB.tr) and re-express the Hamiltonian as 
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Notice that only the symmetric combination of lead op- 
erators couples to the impurity orbital, which plays an 
important role in the solution later. 

In order to work out the flow equation solution for the 
current, it turns out to be convenient to use a finite sys- 
tem with a discrete level spacing. The thermodynamic 
limit will then be taken at the very end when the cur- 
rent is evaluated. We take a constant level spacing A 
corresponding to a constant and equal density of states 
p = 1/A in both leads. The symmetric non-interacting 
terms in the Hamiltonian can then be diagonalized (45j 



by defining the pre- diagonalized basis 



Csa = Z —B s Ck+a + B s d a , 



with the transformation coefficient B s 
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and 



the linewidth T = pirV 2 . The inverse tranformation is 
d a = B s c sr7 and through this the interaction term 
can also be expressed in the pre-diagonalized basis: 



Uru ! n ir = 2J UB sli B Sl B s > 2 B S2 cl,^c Sl ^c 



(5) 
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In the sequel we will work with normal-ordered ex- 
pressions. In this model we define normal ordering with 
respect to the non-interacting ground state in equilib- 
rium, which is also the initial state at time t = 0. The 
chemical potentials of the left and right lead are pl and 
pa, respectively, and V s d = M-R — Ml denotes the voltage 



bias. Strictly speaking, the Fermi function in the pre- 
diagonalized basis has no sharp edge even at zero tem- 
perature due to hybridization. But this effect vanishes in 
the thermodynamic limit and we can use 



n s = {c' sa c sa ) 



with the usual Fermi function 
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In this paper we will generally work at zero temperature 
()3 = oo), the generalization to non-zero temperature is 
st r aight forward . 

The starting point of our analysis is the following 
Hamiltonian 

k sa 
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which corresponds to ([1]) with a single-particle energy 
e d = ~f B 2 s (f L {e s ) + f R {e s )). Notice that the energy 
of the impurity level is then related to the lead chemical 
potentials, i.e. at zero temperature by 



u u v sd , 
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where p = i£i±iLe. The natural parameters in an ex- 
periment are — p, V s( i and U (see Fig. I). For con- 
venience the calculations in this paper will be expressed 
through the parameters pl, Pr and U (or p, V s( i and 
U). However, one can easily solve Eq. [9] to find the cor- 
responding value of p for a given ed — p. Obviously p = 
(or pr = —pL = V s d/2) corresponds to the particle- hole 
symmetric point €d — p = — U/2 (see Fig. I). 

The flow equation approach employs suitable infinites- 
imal unitary transformations in order to diagonalize a 
given many-particle Hamiltonian. Thereby a one param- 
eter family H(B) of unitarily equivalent Hamiltonians is 
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generated, where H(B — 0) is the initial Hamiltonian 
([5]) and H(B — oo) the final diagonal Hamiltonian. Such 
a unitary flow can be generated as the solution of the 
following differential equation 



dH(B) 
dB 



[t}(B),H(B)] 
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where rj(B) is an anti-hermitean operator. Wegner 



showed [38| that the so-called canonical choice rj(B) = 
[H{B),H l - at {B)\, where H lnt (B) the interaction part of 
the Hamiltonian, leads to the required renormalization 
group-like diagonalization scheme. Our key approxima- 
tion will be the restriction to second order in U . In this 
approximation the generator 17(B) — rj^(B) + r/ 2 ) (B) 
takes the following form (for more details see Ref. [45[): 
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where 
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clef 
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The flow of the single-particle energies plays no role in 
the thermodynamic limit if one is interested in impurity 
correlation functions or the current. Therefore the final 
diagonal Hamiltonian takes the following simple form 



H(B 
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Due to spin symmetry we only need to calculate the spin- 
up current I-f since Jf(t) = I±(t). 

The Hamiltonian has been diagonalized by the uni- 
tary transformation U(B) corresponding to the gener- 
ator r)(B) given above. We perform the same unitary 
transformation on the current operator 



dI t (B) 
dB 



[v(B),I t (B)] 



(15) 



Here one should notice that energy-diagonal terms like 

c kt Cs ^ c ki Cs ^ : sti11 
00) but have been neglected in (fl3l) . 
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remain in H(B 
This is permitted since these terms are thermodynami- 
cally irrelevant, that is they vanish in the thermodynamic 
limit. 



III. FLOW OF THE CURRENT OPERATOR 

Clearly, the time evolution generated by (fl"3"|) in the di- 
agonal basis is trivial. However, the price we have to pay 
is to express the observable we are interested in in this 
diagonal basis [44|. Specifically, we look at the current 
operator 1 = 1-^ + 1^, where 



la = {d t N La - d t N R<7 )/2 



iV 



^B„{c\ a Ck-<, - h.c). 
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with the initial condition that I^(B = 0) is given by 
(|14p. In the current operator the anti-symmetric com- 
binations Cfc_| stay invariant under the unitary trans- 
formation, while the commutator of cL. and 77 generates 



higher order terms like 



The commutator 



between this term and r\ feeds back into the coefficient 
of c^. For the lowest order correction with interaction 
(second order in U), the ansatz of the flowing current 
operator looks like 



I t (B) 



^ 7s (B)cl t c fc _ t 

sk 



si Sn safe 



+h.c. 



Substituting this ansatz into Eq. 
lowing flow of parameters: 



(16) 



T5l) one finds the fol- 
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The higher order term in 9bM^|J jS2 is neglected since 
we take only terms up to second order in U into account. 

Next we use the simple time evolution in the diagonal 
basis 

I t (B = oo, t) = e iff (°°)* J t (B = oo)e- iff (°°)* (18) 
leading to 

7*(oo,t) = 7s (oo)e it£ % 



A#f a ' a (oo,t) = Af t s |" 2S2 (oo)e Jt(£s i +£3 2- e3 ^. (19) 
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Next we undo the unitary transformation, that is we in- 
tegrate (fTS"|) from B = 00 with initial conditions (fT9|) to 
5 = 0: 

7 t (0,t) = ^ 7s (0,i)ct t c fc _ t 

s/e 

+ E ^m'C '*) : ^ lt c^c S2 ; : c,_ t 

+h.c. (20) 

Our target is actually 7 S (0, t) in this expression as we will 
find in the next chapter that only J s (0,t) contributes to 
the expectation value of the current. 

The solution of Eq. (fTT]) to the second order in U can 
be written as (see Appendix A) 
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IV. CALCULATION OF THE CURRENT 

At time t = the coupling between the leads and the 
impurity is switched on. The initial state is the non- 
interacting ground state, so the expectation value of the 



current operator can be obtained easily: The quartic 
term in Eq. (|TB| is normal-ordered and does therefore not 
contribute to the expectation value. The time-dependent 
current is expressed as 

7 t (t) = < l t (0,t) > 

ls (Q,t)e-^VB s ^ 
= Re 2^ : ~ (/i(efc) - /fi(efc)). 
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With Eq. (|2"Tj) this gives an explicit expression for the 
current (see Appendix B). The summation over s\ and 
s can be calculated analytically. However, one has to 
be careful since there are poles in the function and the 
summation cannot be simply transformed into a princi- 
pal value integration. We employ the following trick to 
circumvent this problem. For example, when calculating 
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write the expression as 



we introduce a second time t' and 
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Obviously f(t,t) is the original function that we are in- 
terested in and f(t, 0) = 0. Now the pole at e s = D can 
be eliminated by partial differentiation with respect to t'\ 



dt' 
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The poles at e s = can be eliminated likewise (see 
details in Ref. |46|) and the result is ^ s — 

e " fct '- e_rt ' Therefore 
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and the original function follows by integration, / (t, t) = 
Jodt'§{r. The key idea of our method is to introduce the 
additional time parameter t' and to get rid of the poles 
by performing derivatives with respect to t . Afterwards 
one can convert the sum into an integral. Finally one 
perform the integration with respect to t' and gets the 
original function. 



5 



We divide the current into the zeroth order term and 
interaction corrections (see Appendix B), 



I(t) =/(°)(t)+/W(t), 



(27) 



where 
J(°)(t) 

ry/i 



and 
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The dimensionless function T is defined as 



(29) 
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If one uses the hybridization T as the unit of energy and 

1/r as the unit of time, one can write — — as a function 

T/h 

of three dimensionless quantities: V s d = V s d/T, U = U/T 
and t = Tt with 



I 

T/h 



I(i,V sd ,U). 



(31) 



Two limiting cases deserve special attention. First, it 
is straightforward to verify that the current is actually 
zero at t — as required. The calculation of the steady 
state current when t — > oo is also not difficult. The terms 
proportional to e~ Tt vanish in this limit and we find after 
a short calculation: 



I(t) 



2T 



U 2 + r 2 (e 2 + r 2 ) 

2rrU 2 ~, e 2 -T 2 



XJ72 < f dD^l 
e- D 
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V. RELATION BETWEEN THE CURRENT 
AND THE IMPURITY SPECTRAL DENSITY 

Using Green's function methods, the current can be 
expressed by the lesser Green's function as 

J t(*) = -^E Re ( G ^(*'*) - *))> ( 33 ) 

where Cj: a (t,t) = i(d\.(t)cka\{t))a- According to Meir 
and Wingreen |47| . the lesser Green's function is related 
to the retarded impurity Green's function: 



G&(t,t) =j™dt> (gl a (t,t')^=G<(t',t) 

+g< a (t,t')^=G a (t',t) 



(34) 



where 

gi a {t,t>) = -iO{t-t>y^'-V (35) 
9LM = ie^'-^fkc (36) 
are the conduction band Green's functions and 

G<(t,t') = i(4(t')rf t (t)> (37) 
G a (t,t') = i6(t' -t)({d f (t),d\(t')}) . (38) 

are the impurity Green's functions. Eq. (|33[) can there- 
fore be rewritten 

If(t) = J de k (f kL - f kR ) 

xlm / dt'e^-^G^tj) , (39) 



where we have used the relation G a (t,t') = G r *(t',t). 

The retarded Green's function G r (t,t') defined above 
depends not only on the time difference t — t'. We there- 
fore define a time-dependent impurity spectral density 



p(t,e) = — ImG r (t,e), 

TT 

where G r (t, e) is defined via 

POO 

G r (t,e) = / dt'e i < t - t ">G r {t,t') 
Jo 



(40) 



(41) 



Now the time-dependent Meir- Wingreen formula relates 
the time-dependent current with the time-dependent im- 
purity spectral density, 



/(*)= / de(f R (e)-f L (e))p(t,e). 



(42) 



The flow equation result for the Heisenberg time evo- 
lution of dl(t) has already been given in Sect. lIVl There- 
fore the calculation of the time-dependent impurity spec- 
tral density is straightforward, details can be found in 
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Appendix C. Explicit comparison of Eqs. (|28|) and (|29|) 
from the direct solution of the Heisenberg equations of 
motion for the current operator with Eq. (|C6j) shows that 
our previous results in Sect. IIVI are consistent with the 
time-dependent Meir-Wingreen formula as should be ex- 
pected. In the steady state limit t — > oo we find the 
familiar equilibrium impurity spectral density 



lim pit, e) = 



T 2 | 2U 2 eT 
7r(e 2 +T 2 ) + 7r(e 2 +T 2 ) 2 

, i/ 2 f( £ )( £ 2 ~r 2 ) 



dD 



f{D) 



(e 2 + r 2 ) 2 ' 

This equation reproduces the result in Ref. [l5j |. 



D 

(43) 



VI. TIME-DEPENDENT CURRENT AT 
PARTICLE-HOLE SYMMETRY 

The above formulas for time-dependent current and 
spectral density hold for arbitrary left and right lead 
chemical potentials. In the sequel we will present some 
explicit results for the time-dependent current at the 
particle-hole symmetric point, e^ — + fiji)/2 = —U/2. 

We perform numerical integration to get the time- 
dependent current curves. A direct estimation of Eq. ([25)1 
is difficult because there is a pole in the integrand. Alter- 
natively, we calculate the time derivative of the current, 
i.e. 



d_ 

dt 



jW(t) 
T/h 



4U 2 sin %t 



Ft 



Re- 



-iDt 



dD T(D) 

e- rt Tte- rt 
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(44) 
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We then perform numerical integration of the right side 
in (|44|) and employ a fourth-order Runge-Kutta method 
to solve (|44|) and get the current. The symmetry of T 
function, i.e. T(—D) = T(D), is used to simplify the 
calculation. 

Fig. 2 shows the interaction correction to the current 
at different voltage bias. Its time derivative at t = van- 
ishes. This is contrary to the free current, which has a 
sharp increase at t = (see Figs. 3, 4), which indicates 
the initial condition = 0. However, this initial charg- 
ing process is independent of U due to the lack of elec- 
trons in the impurity, which explains 4fl^ c \t = 0) = 0. 

For t ^> 1/r the current correction approaches its 
steady value. Larger voltage bias leads to a stronger 
suppression of the current due to the l/ 2 -dependent cor- 
rection term. This can be understood to arise from 
shot noise decoherence effects, which suppress the quasi- 
particle resonance, similar to the well-established effect 
of current-induced decoherence in the nonequilibrium 
Kondo model 

The suppressed ringing oscillation in both current cor- 
rection and total current can be seen at large voltage bias 
V sd = 2T (see Figs. 3 and 4). From (gSJ and flUJ one 
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FIG. 2: The current correction I^ c \t) due to interaction at 
particle-hole symmetry, = —U/2, for zero temperature. 
The interaction strength is U = T. Results for voltage bias 
Vsd ~ r and Vsd = 2F are depicted. The main features of 
I^ c \t) are a vanishing derivative at t = 0, followed by a sharp 
decrease and finally a smooth crossover towards its steady 
value. One also notices the onset of oscillations at large volt- 
age bias Vsd = 2r. 
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FIG. 3: The current without interaction and for interaction 
strength U = 1.5r at voltage bias Vsd = I\ The interac- 
tion suppresses the current. The inset shows the suppressed 
oscillation of the current. 



can easily deduce the ringing oscillation period 4ir/V sc i, 
consistent with Ref. [48[ . 



VII. CONCLUSIONS 

We have demonstrated how the flow equation method 
(method of infinitesimal unitary transformations) can be 
used to calculate transient and steady state currents in 
and beyond the linear response regime through interact- 
ing quantum impurities. Our approach is perturbativc 
in nature, therefore we are restricted to weak to inter- 
mediate values of the interaction in our analysis of the 
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Integrating with respect to B for a given time t one finds 
the first order solution of M, 



M;lf S2 (B, t) =iVU £ B s[ B 2 S1 B s , B S2 

si 



2(e si + e s < -e Sl -e S2 ) 

ite s , -B(e, +e ,/ -e». -c S9 ) 2 " 

2(6^ +e., 2 -e Sl -e S2 ) 
Taking the limit B — > oo we hnd 
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FIG. 4: The current without interaction and for interaction 
strength U = 1.5F at voltage bias V s d = 2T. The free current 
increases compared to V s d = T, while its interaction suppres- 
sion also becomes stronger due to shot noise induced decoher- 
ence. The inset shows suppressed current oscillation. 
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Anderson impurity model in this paper. One key feature 
of our approach is that there are no secular terms in the 
long time limit, that is the steady state is reached uni- 
formly in the expansion in the interaction. We reproduce 
previous results for the steady state currents [lj| and ob- 
tain analytical results for the transient current behavior 
leading to the steady state. 



Acknowledgments 

We thank M. Moeckel for valuable discussions. We 
acknowledge support through SFB 484 of the Deutsche 
Forschungsgemeinschaft, the Center for NanoScience 
(CeNS) Munich, and the German Excellence Initiative 
via the Nanosystems Initiative Munich (NIM). 



Appendix A: Solution for y s (Q,t) 

The differential equation (| 1 Tf) is solved order by order 
in U. According to the definition of the current opera- 
tor, we have the initial condition 7 S (0,0) = ^Bg and 
M(0,0) = 0. The zeroth order solution can be written 
as M{B, t) = and j s (B,t) = l ^-B s e lUs according to 
Eq. ([T9]). Substituting j s (B,t) into Eq. (Q7]) and inte- 
grating with respect to B at t = 0, we get 



M°£ a \B,0) = iVU ^ 

1 s, s 
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Substituting the above expression and the zeroth order 
solution of 7 S into Eq. (fl~7|) , we find the solution of j s to 
second order in U, 



S ls (t) =7 a (oo,t)-7.(0,t) 
iVBM 2 



si,D 



UD 
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(e s -e sl )(e Sl - D) 



where D — e s j +e s ' 2 — e s 2 an( l T(D) is defined in Eq. 
Then we have 

7.(0, t) = e ^*( 7s (0,0) + <5 7s (0)) - 5 7a (i) 



(A4) 



Si,D 



+ 



{e s - D){e si - D) (e s - e si )(e 51 - D) 

(A5) 



Appendix B: The calculation of the current 

We divide the expression of the current into the ze- 
roth order term and the interaction correction, I^(t) = 

4 0) (t)+4 c) (i), where 



4°>(t) = ReJ2 



2(e s - e fc ) 



(/i(efc) - /fl(efc)), 
(Bl) 



8 



and 



. . iV 2 R 2 TI 2 p~ iekt 



4 c) (t)=Re £ 



s,k,si ,D 



2(e s - e fe ) 



-T(D)B 2 Si 



.(e a -D)(e n -D) (e s - e Sl )(e sl - D) 
X (/i(e*) - /fl(efe)). 

(B2) 

The sum over s and s± is calculated analytically by the 
method introduced in Section IV. The sum over s in the 
zcroth order term is straightforward. Next we need to 
calculate 



A=E 



B: 



-Bi 



efe 



(e s -D)(e sl -D) 



(B3) 



(e s -e si )(e si - Z>) 
We first calculate the sum over s and get 
1 



A =- 



efe - iT 
+ 



gie fc t _ gift e i£ " 



-rt 



D-iY 



-rt 



■zr 



(B4) 



When calculating the sum over s%, we have to get rid 
of the poles at e Sl = D. We rearrange the terms so 
that e Sl — D in the denominator and e lesit — e lDt in the 
numerator appear simultaneously, i.e. 



A : 



1 



i(e 



iDt _ e -rt~) 



2r(L> - ir) 



E B i 



JDt 



(e Sl - D)(e si - efe) 



gift _ gie si t 



(e Sl -D)(e Sl -ir) 



(B5) 



Employing the method from Section IV again we find 



A 



1 



£fc - iT 

+ 



e k - iT e k - D 



D — iT 



g-rt _ £l Dt 



ite 



-rt 



(D - iT) 2 D — iT 



(B6) 



Substituting the expression for A into (|B2j) we obtain an 
expression for I<- c \ The pole at e k = D is a removable 
singularity, so that we can change the sum over k and D 
into a Cauchy principal value integral. This transforma- 
tion makes it easy to estimate the long time limit and to 
compare our result with that in Ref. [15| . The interaction 



correction for the current is then given by 
U 2 f(D) 



I<f\t) = J dDde- 



x Re 



2tt 

ie i(e-D)t _ i 



(/«(e)-Me)) 

(giet-t _ l^ iD + ie _ 2) 



{D-e){D + i) 2 



te l 



(e + i)(D + i) 



(e + i) 2 {D + i) 2 
(B7) 



Appendix C: The calculation of the spectral density 

The evolution of the d\ operator is similar to the cur- 
rent operator and can be expressed as 

4(4) - 5> s (0,4)<4 

S 



&-y a (P,t) and M^ 2S2 (0,4) 



where -%((), 4) 

{{d{t),d\t')}) = J2T s (o,t)Uo,t> 



-^M t |f 2 2 (0, 4). The anticommutator is 



+ E M^ S2 *(0,4)M T ^ 52 (0,4') 

X Qs[s2S' 2 - (C2) 

By using the summation method from the calculation of 
the current, we find 

£ 7s *(0,4) 7s (0,4') 

S 

=e r(t'-t) + u* T ( D ) 

~ 2T(t' - t)e T ^'-^ 
2iT{D + iT) 



e -iDt+iDt' _ e r(t'-t) 



(D + iT) 2 



e —Wt—Tt' _ e iDt'—iDt I e iDt'-Tt _ g - 



r(t+t') 



D 2 + T 2 

(C3) 

Setting B — and performing the summation over s\ in 
Eq. [Ml we get 



(C4) 



AC?? 2 2 (0, 4) = E/-B s ' B s > B S2 — 

& S2 c Sj C S 2 I 6X 

Using the definition D = e s > + e s ^ — e S2 , we obtain 
({d(t),d\t')}) =e r (*'-*) + U 2 T(D) 

' 2T(£_-t)f^^_ 

2iT{D + iT) (C5) 

e -iDt+iDt' _ e r(t'-t) 



(£> + ir) 2 



9 



The impurity orbital spectral density is therefore given 
by 



P(t,e) 



1 



e '(esinei — cosei) 



Tr{e 2 + 1) 



7r(e 2 + 1) 



Re 



U 2 f{D) 



i(e-D)t 



+ 



te 



int-Vt 



{D-e){D + if (D + i)(e + i) 



(e 



iet-Ft 



-l)(iD + ie-2) 



{D + if{e + if 



(C6) 
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